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Abstract 



We have measured the branching fraction and photon energy spectrum for 
the radiative penguin process b — > 57. We find B(b — ► S7) = (3.21 ± 
0.43 ± 0.27±g;i§) x 10 4 , where the errors are statistical, systematic, and 
from theory corrections. We obtain first and second moments of the pho- 
ton energy spectrum above 2.0 GeV, (E 7 ) = 2.346 ± 0.032 ± 0.011 GeV, and 
{Efy - (E y ) 2 = 0.0226±0.0066±0.0020 GeV 2 , where the errors are statistical 
and systematic. From the first moment we obtain (in MS, to order 1/Mg and 
P a 2 s ) the HQET parameter A = 0.35 ± 0.08 ± 0.10 GeV. 
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As a result of several difficult calculations JI|-0 the Standard Model (SM) expression for 
the b — > 57 branching fraction has been determined to next-to-leading order. The numerical 
value originally given [0 was (3.28 ± 0.33) x 10 -4 . Recently, Gambino and Misiak M argue 
for use of a different choice for the charm quark mass, and obtain (3.73 ± 0.30) x 10 -4 . 
The theoretical literature on the connections between b — > 57 and beyond-SM physics is 
extensive [[|, including connections to SUSY, Technicolor, charged Higgs, extra dimensions, 
anomalous WW'y couplings, dark matter, and more. Thus a measurement different from 
the SM value would indicate beyond-SM physics, and a measurement close to the SM value 
would constrain the parameters of SUSY, Technicolor, and other beyond-SM physics options. 

In contrast, the photon energy spectrum is insensitive to beyond-SM Physics. The photon 
mean energy (E^) is, to a good approximation, equal to half the b quark mass m&, while the 
mean square width of the energy distribution, (£?) — (E^) 2 , depends on the mean square 
momentum of the b quark within the B meson. A knowledge of the b quark mass and 
momentum allows a determination of the CKM matrix element V c b from the B semileptonic 
width, as shown in the following Letter ||. Further, the spectrum provides information 
allowing a better determination of the CKM matrix element V u b from the yield of leptons 
near the endpoint of b — > utv |7||||8|]- 

CLEO's previously published measurement || of the b — > 57 branching fraction, (2.32 ± 
0.57±0.35) x 10~ 4 , was based on 3.0 ftT 1 of luminosity, on resonance plus off resonance. The 
photon energy spectrum then obtained was not sufficiently precise to be of interest. Here we 
present a new study of b — > 57, based on 9.1 fb~ on the T(4S) resonance and 4.4 fb~ at a 
center-of-mass energy 60 MeV below the resonance (and below BB production threshold). 
We present the branching fraction and the first and second moments of the photon energy 
spectrum. From the first moment, we obtain |10| , |TTH the HQET parameter A (physically, the 
energy of the light quark and gluon degrees of freedom). 

The characteristic feature of the b — > 57 decay is the high energy photon, roughly mo- 
noenergetic, with E 1 m m^jl ~ 2.3 GeV. In our previous measurement we considered the 
spectrum above 2.2 GeV, incurring a model-dependent correction for the fraction of the 
spectrum below 2.2 GeV. Here we use the spectrum down to 2.0 GeV, which includes ~90% 
of the b — ► sj yield. 

The data used for this analysis were taken with the CLEO detector [12] at the Cornell 
Electron Storage Ring (CESR), a symmetric e + e~ collider. The CLEO detector measures 
charged particles over 95% of 47r steradians with a system of cylindrical drift chambers. 
(For 2/3 of the data used here, the innermost tracking chamber was a 3-layer silicon vertex 
detector Its barrel and endcap Csl electromagnetic calorimeters cover 98% of 4ir. 

The energy resolution near 2.5 GeV in the central angular region, | cos# 7 | < 0.7, is 2.6% 
r.m.s., including a low-side tail. Charged particle species (ti ± , K ± , p/p) are identified by 
specific ionization measurements (dE/dX) in the outermost drift chamber, and by time-of- 
flight counters (ToF) placed just beyond the tracking volume. Muons are identified by their 
ability to penetrate the iron return yoke of the magnet. Electrons are identified by shower 
energy to momentum ratio (E/P), track-cluster matching, dE/dX, and shower shape. 

We select hadronic events containing a high energy photon in the central region of the 
calorimeter (| cos 6^] < 0.7). The photon must not form a 7r° or 77 with any other photon in 
the event. 

There is a large background of high energy photons from continuum processes: initial 
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state radiation; photons from decays of 7r°, rj, and other hadrons. We use several techniques 
to suppress the continuum background, measure what survives with the below-resonance 
data sample, and subtract it from the on-resonance sample. 

There is also a significant background of photons from B decay processes other than 
b — > sj, particularly in the lower portion of our window, 2.0 - 2.2 GeV. We determine the 
dominant background sources, photons from unvetoed 7r° — > 77 and rj — > 77, by directly 
measuring the 7r° and r\ spectra, estimate the K\ and n backgrounds using shower shape, 
and estimate the remaining background sources via Monte Carlo techniques [{14 1 . 



Part of our continuum suppression comes from eight event shape variables: normalized 



Fox- Wolfram |15| second moment R2, S± (a measure of the momentum transverse to the 



photon direction |16|]), R' 2 (the value of R2 in the primed frame, the rest frame of e + e~ 
following an assumed initial state radiation of the high energy photon, with R 2 evaluated 
excluding the photon), cos#' {9' the angle, in the primed frame, between the photon and 
the thrust axis of the rest of the event), and the energies in 20° and 30° cones, parallel and 
antiparallel to the high energy photon direction. While no individual variable has strong 
discrimination power, each posesses some. Consequently, we combine the eight variables 
into a single variable r which tends towards +1 for b — > 57 events and tends towards -1 for 
continuum background events. A neural network is used for this task [|T^] . 

We obtain further continuum suppression from "pseudoreconstruction" , and from the 
presence of leptons. In "pseudoreconstruction" , we search events for combinations of particles 
that reconstruct to a B —>■ X s ^ decay. For X s we use a Kg — > 7r + 7r~ or a charged track 
consistent with a K^, and 1-4 pions, of which at most one may be a 7r°. We calculate the 



candidate B momentum P, energy E, and beam-constrained mass M = J E% eam — P 2 . A 
reconstruction is deemed acceptable if it has x\ < 20, where 



E-E beam \ 2 , ( M — Mb y 2 



xl = [ — I + x 

V o E J \ a M 

where, <je = 40 MeV, % = 4.0 MeV. If an event contains more than one acceptable recon- 
struction, the one with the lowest Xb * s chosen. (It is not important that the reconstruction 
be correct in detail, as we are only using it for continuum suppression. It is correct 40% of 
the time.) 

For events with an acceptable reconstruction, we add Xs> an d | cos# tt |, where 9 tt is the 
angle between the thrust axis of the candidate B and the thrust axis of the rest of the event, 
to the list of variables for continuum suppression. If the event contains a lepton (muon or 
electron), then we add the momentum of the lepton, Pi, and the angle between lepton and 
high energy photon, 9^, to the list. We thus have events with both a pseudoreconstruction 
and a lepton, for which we use r, Xb-> I cos @tt\, Pi, and 9^ for continuum suppression; events 
with only a pseudoreconstruction, for which we use r, x%i an d | cos 9 a \', events with only a 
lepton, for which we use r, Pg, and 6^ 7 ; and events with neither, for which we use r. In each 
of the four cases, the available variables are combined using a neural net. All nets are trained 
on signal and continuum Monte Carlo. We convert the output of each of the four nets, rj, 
into weights, w(rj) = s{rj)/[s{rj) + (1 + a)b{rj)} , where s{rj) is the expected signal yield and 
b(rj) is the expected continuum background yield for that value of the net output r,-, and a 
is the luminosity scale factor between on-resonance and off-resonance data samples (~ 2.0). 
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Weights so defined minimize the statistical error on the b — > 57 yields after off-resonance 
subtraction, as described below. 

Were we to use the weights as just defined, the efficiency would be higher for low mass X s 
states than for high mass X s states. We thus modify the weights for those cases where there 
is a reconstruction, and thus a measurement of Mx a , to render the efficiency independent 
of the X s mass. With these weights, the statistical error is slightly worse, but a systematic 
error has been reduced. The dependence of efficiency on Mx s for events for which there is 
no pseudoreconstruction remains. 

We sum weights on and off resonance, and subtract the off-resonance sum, scaled by 
LjE 2 cm and with momenta scaled by the On/Off energy ratio, from the on- resonance sum 
(the On-Off subtraction). Biases in this subtraction procedure have been estimated with 
Monte Carlo to be (0.5 ± 0.5)% of the subtraction; this has been applied as a correction and 
as a contribution to the systematic error. The on-resonance yields and scaled off-resonance 
yields are shown, as a function of photon energy, in Fig. la, and tabulated in Table 1. The 
On-Off subtracted yields are shown in Fig. lb; they consist of a component from b — > sj, 
and a component from other B decay processes. 

We investigated the component from other B decay processes with a BB Monte Carlo 
sample that included contributions from b —>■ u and b — > sg processes, as well as the dominant 
b — > c decay. We found that the overwhelming source of background (~90%) is photons from 
7T° or 1] decay. Consequently we carefully tuned the Monte Carlo to match the data in 7r° 
and 77 yields, and also rj' and lu. We included radiative ip decays, p — > nj, a x — > ttj, and 
final state radiation. We determine the background from neutral hadrons, K\ and n, using 
information from the measured lateral distribution of the showers. The yield from B decay 
processes is shown in Fig. lb and tabulated in Table 1. The error on this yield is dominated 
by the error on measured On-Off- subtracted n° and rj yields, and by uncertainty in K\ and 
n backgrounds from fits to shower shape. 

The fully subtracted spectrum, On - Off - other B decay processes, is shown in Fig. 2. 
There is a clear peak near 2.3 GeV, indicative of b — > 57. The region of interest for b — *■ S7 
is 2.0 - 2.8 GeV. The region above 3.0 GeV thus serves as a control, indicating how well we 
have subtracted the continuum background, while the region 1.5 - 2.0 GeV serves as another 
control, indicating how well we have accounted for the other B decay processes. The yields 
in both regions are consistent with zero. Subtracting 5% more or less BB background causes 
the B background control region to have a deficit or excess, at the la level. We take ±5% 
of the BB subtraction as the systematic error on the subtraction. 

We determine the efficiency, weights per b — > 57 event generated with photon energy 
above 2.0 GeV, via Monte Carlo simulation. To model the decay b — > sj, i.e., B — > X s ~f, 
we use the spectator model of Ali and Greub [17j , which includes gluon bremsstrahlung and 
higher-order radiative effects. Rather than considering the Fermi momentum parameter Pp 
and the b quark average mass {nib} as quantities given from first principle, we treat them as 
free parameters which allow us to vary the mean and width of the photon energy spectrum. 
We fit our measured spectrum over the range 2.0 - 2.8 GeV, and use the values of the 
parameters so obtained, and their errors, to define our model of B — > X/y. For hadronization 
of sq into X s , we have compared two approaches. In the first, we have taken several K* 
resonances and combined them in proportions to approximate the X s mass distribution 
given by the Ali-Greub model. In the second, we have let JETSET hadronize the sq, again 
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giving the Ali-Greub X s mass distribution. We include in the systematic errors on the 
efficiency an uncertainty in the modelling of the B — > X s 7 decay (Pf, (m>b), hadronization) , 
an uncertainty in the simulation of the detector performance (track-finding, photon-finding, 
resolutions), and an uncertainty in the modelling of the other B. 

As an alternative to the Ali-Greub model for B — > X s 7 decay, we have used the descrip- 
tion of Kagan and Neubert ||, with their simple two-parameter shape function, with the 
two parameters m b and a adjusted to fit our measured photon spectrum. We obtain results 
very similar to those obtained using Ali-Greub, both for efficiency and for the moments. 
(Because we determine our B — > X s 7 model from a fit to data, our results are independent 
of the nominal description, as long as that description allows us two adjustable parameters 
to match data on mean and width.) 

To obtain the b —>■ sj branching fraction, we take the yield between 2.0 and 2.7 GeV, 233.6 
± 31.2 ± 13.4 weights, where the first error is statistical and the second is systematic, from 
BB (± 5%) and continuum (± 0.5%) subtraction. The efficiency is (3.93±0.15±0.17) x 10~ 2 
weights per event, where the first error is from model dependence of the B — > X s 7 decay and 
the second error is from detector simulation and model dependence of the decay of the other 
B. Our sample contains 9.70 million BB pairs (±2%). We obtain an uncorrected branching 
fraction of (3.06 ± 0.41 ± 0.26) x 10" 4 . 

This uncorrected branching fraction, what we directly measure, is the branching fraction 
for b — * 57 plus b — > dj, for 5-rest-frame photon energies above 2.0 GeV. We apply two 
"theory" corrections to obtain a branching fraction of more direct interest. Using a model for 
b —> d*y similar to that used for b — > 57, we find the efficiency for b — > d'-f to be the same as 
that for b — > S7. (Pseudoreconstruction favors b — > 57 over b — > d^ but shape variables and 
photon energy spectrum favor b — > d'y over b — > 57.) The SM expectation is that the b — *> d'y 
and b — > S7 branching fractions are in the ratio |V^/Vi S | 2 . Using IV^/V^I = 0.20 ± 0.04 |18| 
(taking V ts = - V c b, from unitarity), we correct the branching fraction down by (4.0 ± 1.6)% 



of itself, to remove the b — > c/7 contribution. The fraction of b — > 37 decays with photon 
energies above 2.0 GeV is sensitive to the b mass and Fermi momentum. We use 0.915to;o55, 
as given by Kagan and Neubert [p|Jl9l to extrapolate our branching fraction to the full energy 
range (actually, to energies above 0.25 GeV). With these two corrections, we have 



B(b -► s 7 ) = (3.21 ± 0.43 ± 0.27±^) x 10" 4 



for the branching fraction for b — > S7 alone, over all energies. The first error is statistical, the 
second systematic, the third from the theory corrections. This result is in good agreement 
with the Standard Model prediction. It is also in agreement with our previously-measured 
result 0, which it supercedes. 

We have calculated the first and second moments of the photon energy spectrum, in two 
ways. In the first, we correct the raw spectrum for the energy dependence of the efficiency, 
calculate moments from that spectrum, apply corrections for experimental resolution and 
for the transformation from lab frame to B rest frame. We then apply a final, empirical 
correction, obtained by carrying out the above procedure on Monte Carlo. In the second 
way, we take our best fit Monte Carlo model, with Pp and (mb) determined from a fit 
to the spectrum, and take the moments given by that model. The two ways of obtaining 
moments agree, and we take their average. We thus obtain moments in the B rest frame, 
for E^ir est frame) > 2.0 GeV: 
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2.346 ±0.032 ±0.011 GeV . 

= 0.0226 ±0.0066 ±0.0020 GeV 2 



Heavy Quark Effective Theory and the Operator Product Expansion, with the assumption 
of parton-hadron duality, allows inclusive observables to be written as double expansions in 



powers of a s and 1/Mb- The expressions [1C], [11] for the moments of the photon energy 
spectrum in B — ► X s 7, for > 2.0 GeV, in the MS renormalization scheme, to order /3o« 2 
and 1/Mj,, are given in Eqs. [1] and ^[ 



(E y ) = - -385f - 

13pi-33p2 



.620/5, 



12A/|~ 



01 



Ti+3T 2 +r 3 - 



-37i P2C2 



9M B Mf,C 7 



- 1.175A)( 
C(l/M|)], 



((JE, - (£ 7 » 2 > = M\ [5=^ + (0.00815^ + 0.01024/5 (f ) 2 ) 
- ^(0.05083^ + 0.05412/3 (t) 2 ) - 2 §Sf - l£f + 0{l/M%)i (2) 

where C2 and Cj are Wilson coefficients and /3q is the one-loop QCD (3 function. The 
parameters pj, are estimated, from dimensional considerations, to be ~ (0.5GeV) 3 . Using 
Eq. |l], we obtain 

A = 0.35 ±0.08 ±0.10 GeV , 

where the first error is from the experimental error in the determination of the first moment, 
and the second error is from the theoretical expression, in particular from the neglected 
1/Mg terms, and the uncertainty of the scale to use for a s , which we take to be from m&/2 
to 2m&. 

The expression for the second moment converges slowly in 1/Mb, and so we make no 
attempt to extract expansion parameters from it. 

In summary, we have measured the branching fraction for b — ► 57 to be (3.21±0.53) x 10~ 4 , 
in good agreement with the Standard Model expectation. We have measured the first and 
second moments of the B rest frame photon energy spectrum above 2.0 GeV to be (E^) = 
2.346 ± 0.034 GeV , and (£ 2 ) - (£ 7 ) 2 = 0.0226 ± 0.0069 GeV 2 . From the first moment, 
we have obtained a value for the HQET OPE parameter A = 0.35 ± 0.08 ± 0.10 GeV . 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. We thank A. Falk, A. Kagan and M. Neubert for many 
informative conversations and correspondences. This work was supported by the National 
Science Foundation, the U.S. Department of Energy, the Research Corporation, and the 
Texas Advanced Research Program. 
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Energy Range 


1.5 - 2.0 GeV 


2.0 - 2.7 GeV 


3.0 - 5.0 GeV 


On 


2718.8± 17.2 


1861. 7± 16.5 


1083.2± 8.0 


oxOff 


1665. 6± 14.6 


1399.4± 15.31 


1101.4± 11.6 


On-axOff 


1053.2± 22.6 


462. 3± 22.5 


-18.2± 14.1 


B backgrounds 


1033.0± 35.4 


228.6± 21.6 


1.5 ± 1.1 


On -oxOff -B 


20.2±42.0 


233.6±31.2 


-19.8± 14.2 



TABLE I. Yields (summed weights) with statistical errors for three photon energy intervals. 
Given are yields On T(45) resonance, scaled Off-resonance yields, On minus scaled Off, estimated 
backgrounds from B decay processes other than b — > sj and b — > dj, and On minus scaled Off 
minus B backgrounds, the putative b — * 57 plus b — > dj signal. 



1850801*06 




Ey (GeV) 

FIG. 1. Photon energy spectra (weights per 100 MeV). The upper plot (a) shows the On T(4S) 
and scaled Off-resonance spectra. The lower plot (b) shows their difference, and the spectrum 
estimated for B decay processes other than b — > 57 and b — > dj. 
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\ Data 

— Spectator Model 




FIG. 2. Observed laboratory frame photon energy spectrum (weights per 100 MeV) for On 
minus scaled Off minus B backgrounds, the putative b — > S7 plus b — ► signal. No corrections have 
been applied for resolution or efficiency. Also shown is the spectrum from Monte Carlo simulation 
of the Ali-Greub spectator model with parameters (mi,) = 4.690 GeV, Pp = 410 MeV/c, a good 
fit to the data. 
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